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ABSTRACT 

\Q With respect to the two theoretical approaches for estimating the hadronic 

ON productions of the mesons B c and B*: the full a 4 calculation and the 

fragmentation approximation, with a thorough comparative study a quite 
remarkable conclusion may be drawn. In hadronic collisions the condition for the 
Qh applicability of the fragmentation approximation is the transverse momentum 

Pt of the produced B c and B* being much greater than the meson masses, and 
the higher the energy of the relevant subprocess the higher the Pt is requested. 
The mechanisms for the productions are also clarified. In order to show the 
resultant differences of the two approaches, newly introduced distribution of the 
fraction of the B c (or B*) meson energy (a measurable observable) in the center 
£2 °f mass of the subprocess is emphasized. 

QN The reason why we would like to talk about the subject here is three-fold: first 

S of all, the problem is interesting itself; secondly, there were some obscurities due to 

G\ a few incorrect numerical results and misleadings, and now the situation has been 

clarified; the third, some fresh results have been achieved. 

The heavy flavored be meson states have attracted considerable interest due to 
^ their interesting properties ( the spectroscopy of be meson states 1 ; the weak decays 2 

4^ and the productions 3-18 etc ). The experimental search for these mesons is now under 

way at high energy colliders such as the LEP e + e~ collider (at Z° resonance) and the 
Tevatron pp collider (at the full energy 1.8 TeV). Various experimental results are 
expected to come out soon thus to estimate the production cross sections as precisely 
as possible has become a desirous theoretical task. Fortunately the mesons are weakly 
bound states and relatively simple, their hadronization is calculable in the framework 
of perturbative QCD (PQCD) in terms of the wave function obtained by potential 
model. At LEP being relatively simple, the productions of the pseudoscalar ground 
state B c and the vector meson state B* are dominated by Z° decay into a bb pair, 
followed by the 'fragmentation' of a b quark into the B c or B* c meson 3 ' 4 ' 5 , whereas 
as first pointed out in Ref.[6], the energetic hadronic productions are much more 
complicated even only to the lowerst order of PQCD in a 4 and dominated by the 
subprocess of 'gluon-gluon fusion' (gg — > B c (B*)bc). In practice, a possible and 
alternative way to calculate the hadronic productions is to apply the fragmentation 
approximation. With the approximation the calculation can then be considerably 
simplified, as indicated in Ref.[7]. Since then one interesting question is addressed, 
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how well and/or in which kinematic region for the hadronic productions of the B c and 
B* mesons the fragmentation approximation works. Therefore several groups, based 
on the same consideration, have recalculated the productions. Namely since Ref.[6] 
presented the numerical results for the hadronic productions first, recalculations have 
been completed and distributed 8-15 . The calculations for the hadronic productions 
of the mesons B c and B* c to order a\ in PQCD involve very complicated numerical 
calculations. At the first stage, not all the calculations were in agreement: in Ref.[8] 
an order of magnitude larger result than Ref.[6] is claimed; in Ref.[9] again a result 
larger than Ref.[6] but smaller than Ref.[8] is obtained; a result similar to Ref.[9] is 
found by Masetti and Sartogo 10 . But more recently, it is found by the same authors 

of Ref.[9] that a color factor - had been overlooked in their previous work and after 
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correcting this factor, a result in agreement with Ref.[6] has been achieved 11 . In 
Ref.[12] similar numerical results to Ref.[6] have also independently been obtained. 
In Ref.[13] different numerical results have been presented, whereas it is difficult 
to compare directly with others due to different parton distribution functions and 
different energy scale being adopted in the calculation. Thus in Refs. [14,15], the 
authors of Ref.[6] as well as their collaborators have re-studied the productions 
extensively with the updated parton distribution functions 19 and made thorough 
comparisons with others' 11_13 . Now at least of the four groups 6 ' 14 ' 15 ' 11 ' 12 ' 13 , the 
numerical results for the cross sections of the subprocess are consistent each others 
within the uncertainties of the numerical calculations. One may now be confident 
that the numerical results of the full order a\ PQCD calculations are in agreement. 
Thus I will discuss the problem based on the confident results. 

From experience, it is naively believed that, when the transverse momentum 
P T of the meson B c or B* is large, the hadronic productions of the mesons are 
dominated by jet fragmentation, so the fragmentation approximation may become 
very valid. However one should carefully examine the approximation not only because 
it will potentially be used as a further test of the factorization theorem of QCD but 
also because of special interest in estimating the productions as reliable as possible 
i.e. as mentioned above to predict the productions precisely so as to guide the 
discovery of the meson B c in experiments. To pursue the goal, subsequent comparisons 
between the PQCD full a 4 calculation and the fragmentation approximation were 
made by several groups 6 ' 14 ' 15 ' 11 ' 13 . Though all the results were comparable, different 
conclusions still were drawn due to different observables in different kenetics regions 
being emphasized. I will talk about the aspects in detail below. 

Based on PQCD to calculate the hadronic productions, the full a A s order calculation 
is illustrated in Fig. la with Fig. lb and the fragmentation approach in Fig. la and 
Fig.lc. 

As pointed out above, it is very interesting to investigate the appoaches. Since 
the P T distribution decreases very rapidly as P T increases, one would lose a lot of 
statistics if one had only considered large Pt events. Thus the comparitively low 
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P T events should also be considered carefully, as long as they may escape from the 
experimentally proper cuts. In Ref.[ll] by examining the production ratio for B* to 
B c , it is claimed that the fragmentation approximation breaks down even for very 
large Pt, whereas in Ref.[13] by investigating the Pt distribution of the B c meson, 
it is claimed that the fragmentation approximation works well if Pt exceeds about 
5 — 10 GeV. Whereas we investigated this problem quite early 14 and re-examined it 
more carefully thus a more completed feature has been achieved 15 . 




Figure 1: The factorizations for the productions 

To clarify the problem, let us focus the discussions on the subprocess. To the 
lowest order (a 4 ), there are 36 Feynman diagrams responsible for the dominant gluon 
fusion subprocess g(ki) + g(k 2 ) — >■ B c (p) + b(q 2 ) + c(qi), where hi, k 2 , p, qi, and q 2 are 
the respective momenta. As pointed out in Refs.[6,14], of the 36 Feynman diagrams 
one may split them into five independent "groups" according to the color structure 
and each of them is gauge invariance for QCD. Furthermore, the contributions from 
the kinematic region where certain factors of the amplitude for the process are nearly 
singular; i.e., some of the internal quark lines and gluon lines in the Feynman diagrams 
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are close to mass-shell are substantial even dominant, especially when the c.m.s. 
energy of the subprocess, Vs, is much larger than the heavy quark mass. Specifically, 
here for the concerned subprocess and in a special chosen gauge where large number 
cancellation due to gauge invariance does not happen at all, the possible singularities 
may arise from the inverse power(s) of the following factors, or their products: 

qi-kj, p-kj, + and (kj - a x p - q^ 2 , (1) 

which corresponds to the denominators of the quark and gluon propagators appearing 

hi t i t 

in the amplitude and here i,j = 1,2 and a\ 2 = 7 r is the ratio of quark masses. 

(m c + m h ) 

The singularities result in the cross section, for the subprocess and upto the lowest 

I p 

twist contributions, proportional to - 2 C (where fs c is the B c decay constant) and 

s M Bc 

with some logarithmic correction terms such as \n(s/M Bc ) being involved. When the 
P T of the B c meson is large, only (oyj + qi) 2 can still be small (~ m 2 ), therefore the 
fragmentation functions can then be extracted from the most singular part containing 
the inverse powers of this factor in the square of the amplitude. It then follows that in 
the large Pt region the subprocess is dominated by the fragmentation. However, when 
the Pt of the B c meson is small, the produced B c , as well as the b and the c quarks, 
can be soft or collinear with the beam. In this region the amplitude is highly singular 
because two or more of the internal quarks or gluons in certain Feynman diagrams can 
simultaneously be nearly on-mass-shell. Although this region is a smaller part of the 
phase space, these nearly singular Feynman diagrams make a substantial contribution 
and dominate to the cross section in the region. In fact, in the square of the amplitude 
we can isolate all the terms which contribute to the lowest twist cross sections using 
singularity power counting rules 20 . When \/I ^> M Bc the lowest twist contributions 
dominate, while the higher twist contributions are suppressed by a factor m 2 /s at 
least and small. We can decompose the terms which contribute to the lowest twist 
cross sections into three components: One is due to the fragmentation contribution, 
which dominates the large Pt region, and the second is due to the non-fragmentation, 
which comes from the other singular parts as discussed above (those in which two or 
more quarks or gluons are nearly on-mass-shell) . The non-fragmentation dominates in 
the smaller P T region. As for the third, the regular parts, they contribute the smooth 
"background" only and small, thus they are not so interesting. The contributions 
from the first two components are quite clearly distinguishable in the Pt distribution 
of the subprocess, particularly at large Vs. 

In order to show the factor quantatively, the cross sections for the dominant 
subprocess gg — >■ B C (B* C ) + • • • versus P T (without any P T cut) are calculated by 
the two approaches and at four energies for the colliding gluons: 30GeV, 60GeV, 
lOOGeV and 200GeV. The results are plotted in Fig. 2 (Fig. 2a-2d). Here the solid 
lines indicate those for the full QCD a A s calculations; the dotted lines indicate those 
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for the fragmentation approximation. Note that in the calculations here the values 
of the parameters a s = 0.2, mj = 4.9 GeV, m c = 1.5 GeV, M Bc = 6.4 GeV and 
Jb c = 480 MeV are taken, moreover in the fragmentation calculations, in order to 
reduce the error caused by the phase space integrations, we directly used the squared 
matrix elements from which the fragmentation functions are derived, rather than the 
fragmentation functions themselves. 




O lO 20 30 +0 50 30 40 SO SO 1O0 

P T (GeV) P T (GeV) 

Fig Bc Fig 3d 



Figure 2: The cross sections for various colliding energies 

Let us take the result at a very high colliding enery = 200 GeV as an 
example to compare the full a* calculation and the fragmentation approximation. 
It is easy to see from Fig. 2d that when P T is larger than about 30 GeV for 
the B c and about 40 GeV for the £?*, the fragmentation approximation is close 
to the full a\ calculation. However, when the value of the P T becomes smaller 
and smaller than about 30 GeV for the B c and about 40 GeV for the B*, the 
deviation between the fragmentation approximation and the full calculation becomes 
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larger and larger. The non-fragmentation component clearly starts to dominate the 
productions at certain critical value of Pt, which certainly is much larger than the 
B c meson mass. If comparing all the four figures Fig. 2a-2d it is easy to see that 
the critical value of Pt increases slowly with s is increasing, that indicates there is 
an additional enhancement maybe due to logarithmic terms such as lns/m 2 in the 
non- fragmentation component compared to the fragmentation component. When 
is not very large this two component decomposition is less distinct, since the higher 
twist terms cannot be ignored. In summary, the fragmentation is not always to be a 
very good approximation 14 ' 15 ' 11 ' 13 . 

It is worth while here to talk a little about the interesting and similar process, the 
production of B c and B* mesons in photon-photon collisions. A comparative study of 
the full o?o? s calculation with the fragmentation approximation in the photon-photon 
process was presented in Ref.[16,17], where it was claimed that the fragmentation 
approximation is not valid. There are 20 Feynman diagrams which can be divided 
into four gauge invariant subsets corresponding to various attachments of photons 
onto the quark lines; i.e., subsets I, II, III, and IV corresponding respectively to the 
attachment of both photons onto the b quark line, onto the c quark line, one photon 
onto the b quark line and the other onto the c, and the interchange of b and c. Subset 

I is dominated by the b quark fragmentation into the B c when P T is large, as discussed 
above. Subsets III and IV, the so-called recombination diagrams, can contribute to the 
non-fragmentation component substantially and decrease rapidly as Pt is increasing, 
especially, when P T ^> M Bc . However, subset II is somewhat unusual. Although 
this contribution is relatively suppressed by the smaller probability for subsequent 
bb quark creation, it nevertheless gives quite a large contribution to the total cross 
section because of the enhancement of the c quark electric charge. For instance, it 
has been found that when = 100 GeV the result of the full calculation is an order 
of magnitude larger than the fragmentation calculation, even for large Pt- However, 
when \/I becomes extremely large; e.g., = 800 GeV, the contribution of this subset 

II is dominated by the c quark fragmentating into the B c meson when P T ^> M Bc . 
This implies that naive power counting rules are violated 20 . However, such a kind of 
contribution in photon-photon collisions, which involves b, b quark pair creation, is 
not important in the hadronic productions of the B c and B*, simply because such an 
enhancement factor of 16, the ratio of the electric charges of c and b quarks, does not 
occur at all. 

The factors emphasized above may also be shown in the way: to present the 
cross sections for the hadronic productions of the B c and the B* mesons with various 
Pt cuts, that is related to experiments more directly. Considering the situation of 
the experiments and experimental results at Tevatron being available soon, the cross 
sections at Tevatron energy and with various cuts for P T and a fixed one |F| < 1.5 
are calculated by the two approaches: the a 4 calculation and the fragmentation 
approximation. The results are shown in Table I. Here the CTEQ3M parton distribution 
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functions 19 is adopted. 



Table I. Total 


cross 


sections 


°( P TB C 


> PTmin) at Tevatron in 


nb 




PTmin (GeV) 





5 


10 


15 


20 


25 


30 




°B e {<4) 


1.8 


0.57 


0.087 


0.018 


4.8 x 10~ 3 


1.6 x 10~ 3 


6.3 x 10" 


-4 


VB c (frag.) 


1.4 


0.47 


0.071 


0.014 


4.0 x 10" 3 


1.3 x 10" 3 


5.3 x 10" 


-4 


a B * c (a A s ) 


4.4 


1.4 


0.22 


0.041 


1.1 x 10~ 2 


3.4 x 10~ 3 


1.3 x 10" 


-3 


°b* c ( frag. ) 


2.3 


0.78 


0.12 


0.025 


6.8 x 10~ 3 


2.3 x 10~ 3 


9.2 x 10" 


-4 


cb* ( frag. ) 


0.52 


0.55 


0.56 


0.61 


0.63 


0.67 


0.70 





We should note here that the values in the table are greater than those obtained 
in the Refs.[6,14] by a factor about two, it is due to more careful and reasonable 
considerations on the coupling constant a s i.e. in the estimates for Table I, the 
coupling constant a s (P^) (with P T ^> mj) is adopted when it related to the part of 
the b and b pair production (with a high P T in the considering mechanism), whereas 
a s {Am 2 c ) is adopted for the part of the fragmentation (the cc pair producion and the 
hadronization for the meson). 

As pointed out by us 14 ' 15 , it is difficult to make a definite conclusion only from 
the Pt distribution consideration if and where the fragmentation approximation is 
reliable for the hadronic productions of the B c and B* mesons. In fact, we have found 
that sometimes the consideration of P T alone can be misleading. 

To clarify this issue, the newly introduced distribution a(z) 15 : 



a(z) = J dxidx 2 gi{x u n)g 2 {x 2 , n) 
2{h + k 2 ) -p 



da(V§, /i) 



dz 



(2) 



is much more efficient. Here z = v 1 — ^ — — and gAxi, /i) are the gluon distribution 

s 

functions. In the subprocess center of mass z is simply twice the fraction of the total 
energy carried by the B c or B* meson. The distribution a(z) provides a sensitive 
means to investigate the dynamics of the production processes and the fragmentation 

approximation. Clearly, if the fragmentation approximation is valid, ^ 
be factorized further as 



dz 



can 



da(y / li, n) 



dz 



(3) 



where Dq^bXz, A*) are the usual fragmentation functions and ^ gg ^Q i Q i is the gluon 
fusion subprocess cross section for the production of the heavy quark pair QiQi. In the 
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approximation, the integrals over x\ and xi can be performed, and the fragmentation 
function can be factored out: a(z) is simply related to the fragmentation functions 
but not much sensitive to the parton distribution functions and the kinematic cuts as 
well. Therefore, comparing a(z) calculated in the fragmentation approximation with 
the full order a\ calculation, provides a quantitative criterion to judge the validity 
of the fragmentation approximation. Note that z is a very useful variable and is an 
experimentally measurable quantity too. 





Figure 3: da/dz versus z at Tevatron 



In Fig 3, the two figures: Fig 3a with a cut Pt > 10 GeV and Fig 3b with cuts 
Pt > 20 GeV and P T > 30 GeV are included, and the solid lines indicate those for 
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the full QCD at A s calculations; the dotted lines indicate those for the fragmentation 
calculations. Since the cross sections for the B c and B* c meson productions decrease 
very rapidly as Pt increases, the distribution a(z) is sensitive to the smallest Pt 
region through given P T cuts. From Fig. 3, some general features are evident: as for 
the B c meson, the distribution a(z) obtained by the fragmentation approximation is 
overestimated in the higher z region while it is underestimated in the lower z region 
for a small Pt cut; but, after integration over z, the result is similar to the full a\ 
calculation, whereas, as for the B* meson, even for a large P T cut, the distribution a(z) 
is underestimated at all values of z and, after integration over z, the result is definitely 
smaller than the full a\ calculation. This shows that it is simply fortuitous that the 
Pt distribution of the B c calculated in the fragmentation approximation is similar to 
that from the full ol\ calculation for Pt below a certain value, particularly down to 
p T M Bc ■ It is also clear that when P T is increased the distributions obtained by the 
two approaches become closer. As shown in Fig. 3b, when P T becomes as large as 30 
GeV the curves calculated in the fragmentation approximation are quite close to the 
full calculation. This indicates that the fragmentation approximation is valid in 
the large P T region and the higher the energy of the subprocess the higher the P T is 
requested. Perhaps the complicated feature pointed out above causes the authors of 
Refs.[ll,13] to draw different conclusions. We should emphasize here one more point 
that the difference between the fragmentation approximation and full calculation is 
not universal, but is process-dependent, and it does not satisfy the Atarelli-Parisi 
evolution. 

Finally, as discussed above, for B* meson production the fragmentation approach 
always underestimates the full a*. The deviation from the full calculation for the B* c 
meson can be used as a criterion to test the validity of the fragmentation approximation 
The results for the total cross section g(Pt > Pt mm) f° r various Pt cuts are listed 
in Table I. Taking agreement within 30% as the criterion for the validity of the 
fragmentation approximation from Table I, one may see that Pt should exceed about 
30 GeV at Tevatron, a value considerably much larger than the heavy quark masses. 
If comparing the ratio of B* to B c production predicted by the full a\ calculation 
with the fragmentation approximation, a similar conclusion may also be reached. We 
note that this conclusion is also rather insensitive to the choice of the QCD energy 
scale fi and the parton distribution functions. 
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